Dispersion of a single hole in the t — J model 
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The dispersion of a single hole in the t — J model obtained by the exact result of 32 sites and 
the results obtained by self-consistent Born approximation and the Green function Monte Carlo 
method can be simply derived by a mean-field theory with d-RVB and antiferromagnetic order 
parameters. In addition, it offers a simple explanation for the difference observed between those 
results. The presence of the extended van Hove region at (tt, 0) is a consequence of the d-RVB 
pairing independent of the antiferromagnetic order. Results including t' and t" are also presented 
and explained consistently in a similar way. 



PACS numbers; 74.20.-z, 71.27.-|-a, 74.25.Fy 

Important advances in the study of high tempera- 
ture superconductors have been made recently by angle- 
resolved photoemission (ARPES) experiments Ex- 
tended van Hove singularity (EVHS) near the Fermi 
surface of the superconductors is found. In particular, 
Wells et al. Q have measured the ARPES for the insu- 
lating Sr2Cu02Cl2 and find good agreement with the 
prediction of the t — J model about the bandwidth and 
other features. Although there is substantial disagree- 
ment about the position of the energy level near (7r,0). 

There were many theoretical studies [Q of the proper- 
ties of a single hole in a quantum antiferromagnet. The 
theoretical data used in comparison with ARPES data 
are obtained from the t — J model by using the self- 
consistent Born approximation (SCBA) to treat the 
scattering of a single hole with the spin waves of the 
Neel state. The result of SCBA shows an EVHS near 
(7r,0). The SCBA dispersion of the single hole is also 
reproduced in the Green function Monte Carlo (GFMC) 
approach by Dagotto et al. Based on these results, 
Dagotto and collaborators ||] have argued that the EVHS 
is due to the antiferromagnetism and have proposed the 
antiferromagnetic-van Hove theory to explain the mech- 
anism of superconductivity. 

Recently Leung and Gooding have solved the spec- 
tra exactly for a single hole in a 32-site square cluster. 
The dispersion relation obtained is very similar to the 
result of SCBA except they have not observed the exact 
symmetry resulted from folding the Brillouin zone (BZ) 
into half due to the antiferromagnetic long range order. 
It should be noticed that this asymmetry is also observed 
in the GFMC result § for a 12 x 12 cluster. 

In this paper we will show that the results obtained 
by SCBA and GFMC [D as well as exact diagonaliza- 
tion [ pO{ can be easily reproduced by a mean-field theory 
including both d-wave resonating valence bond (d-RVB) 
and antiferromagnetic order parameters. The presence 
of the EVHS is associated with the d-RVB pairing. 

We have also looked at a more general t — J model by 



including next-nearest-neighbor hopping, t' , and next- 
next-nearest-neighbor hopping, t" . It is shown below 
that the complicated dispersions obtained by including 
t' and t" are easily understood in terms of the mean-field 
theory. A better agreement with the ARPES of the insu- 



lating Sr2Cu02Cl2 
values of t' and t" ill 



could be obtained by tunning the 
^f2t, but there is still discrepancy. 
Several years ago several groups ||l^-|l^ have used the 
variational Monte Carlo method to study the phase dia- 
gram of the t-J model at small hole density. The antifer- 
romagnetic phase boundary obtained agrees fairly well 
with experiments Near half-filling the best trial 

function has antiferromagnetism and d-wave resonating- 
valence-bond (d-RVB) |l^] coexist. But no analytical 
mean-field theory was discussed. 

The two order parameters, Neel order and singlet 
d-RVB pairing, can be considered simultaneously in a 
mean-field theory of the t — J model, which is given by 
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where {i, j) is the nearest neighbor pairs and Ci^- — Ci^{l — 
ni _fj). At half-filling, we shall only consider the Sj ■ Sj 
term which involves three mean-field order parameters. 
The staggered magnetization is to^ = {S^) = —{Sg), 
where the lattice is divided into A and B sublattices. 
The uniform bond order parameters is x = (So- ^tcr^3<^)^ 
and d-RVB is A = (cjiQi — Cj^Cii) * ^'^^ j ^'^^ nearest 
neighbor sites in x-direction and —A for y-direction. The 
mean-field Hamiltonian can be diagonalized and this is 
done recently by Inaba et al. jl^ in the slave-boson for- 
malism. They have examined the mean-field phase dia- 
gram of these order parameters. We shall adopt a slightly 
different approach. 

Instead of taking into account all three order param- 
eters together in the mean-field Hamiltonian, we first 
consider the staggered magnetization rus and uniform 
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bond order x- They produce upper and lower spin- 
density-wave (SDW) bands with dispersions: ±^fc = 
±(e| + {JmsY)^ where = —jJx{coskx + cosky) The 
states in these bands are then paired. At half filling the 
mean-field Hamiltonian, in addition to a constant, is of 
the form 

Hmf = ^ -Cfeflk^ak^T + X! ^fe^kCT^kCT 

k,cj k,fT 

+ ^Afe(a_k|ak| - fe_k|6kj^) + /i.e. (2) 

k 

where = | JAdfc, and = cosk^ — cosky. The sum is 
taken over the sublattice BZ (SBZ). The operators of the 
lower and upper SDW bands are related to the original 
c operators by: Okcr = a^Ckcr + o'/3kCk-i-Qo-, and b^^ = 
—aPkCka + QfkCk-i-Qo- respectively. We set Q = (tt, tt) 
for the commensurate SDW state, = ^(1 — |^), and 

Hmf can be diagonalized separately for the lower and 
upper SDW bands. For the lower band, in addition to 
a constant, it becomes — £^fe(/ik./ik — /2k /2k) where 
/ik = Wfeaki - Vka\^, /2k = WfcflkT + Uka\^, and 

Ek - (el + Al)i = (el + {Jm^f + Alf^ . (3) 

The coherence factors are = 5(1 + bt)' ^^'^ ^fc ~ 
^(1 -f -^). In terms of the operators aj^, the ground 
state has the familiar ECS form. Similarly, the upper 
SDW band also forms two bands with identical dispersion 
±Ek as those of /ik and /2k- If we had not chosen zero 
chemical potential for the half-filled case, there would be 
four non-degenerate bands as shown by Inaba et al. . 

In the above approach, the local constraint of no dou- 
bly occupied sites is satisfied for the whole lattice on the 
average but it is not for each individual site. To obtain a 
more accurate quantitative result for the Hamiltonian H, 
we shall use the variational Monte Carlo (VMC) method 
that satisfies the constraint exactly. The ground state 
of Hmf is simply given by the product of two BCS-like 
wave functions, one for the lower SDW band and the 
other for the upper SDW band. In the presence of the 
constraint the wave function has the form 

1*0) = Pd{J2{A,a+^a\^ + B,b+^b\^)f^/'\0) (4) 

k 

where is the total number of electrons and coefficients 
Ak = {Ek + a)/Afc and Bk = -{Ek - CO/Afc. The pro- 
jection operator Pd enforces the constraint of no doubly 
occupied sites. In this wave function there are two vari- 
ational parameters: A/x and rus/x- the absence of 
staggered order, m^, this is exactly the same RVB wave 
function used by Gros [Q. Without pairing this wave 
function describes the SDW state. |\I'o) is similar to the 
trial wave function used by Chen et al. but with a 



slightly lower energy. It has about — 0.332J per bond 
which is within one percent of the best estimate of the 
ground state energy of the Heisenberg model. The suc- 
cess of l^'o) gives support to the mean- field theory that 
derives Eq. (4). 

According to the mean-field Hamiltonian Hmf dis- 
cussed above, creating a hole is to take away a quasi- 
particle /ik from the lower SDW band in the ground 
state, or the corresponding one in the upper SDW band. 
Hence the energy of such a state is just —Eg + Et, where 
—Eg is the ground state energy at half-filling. Interest- 
ingly, the quasiparticle dispersion actually has the simi- 
lar form as the dispersion of the hole obtained by SCBA. 
In Fig.l we show that the energy dispersion of a single 
hole obtained by Liu and Manousakis by using SCBA 
for a 20 X 20 lattice with J = 0.2 can be fitted quite 
well by Ek — Eq with parameters, x = 5.47, rus = 13.2 
and A = 2.27. Eq = 5.37 is just a constant shift. In 
this paper, the energy unit is t. However, the value of 
rus is unphysically large. This is mainly due to our ne- 
glect of the constraint in deriving Eq. (3). The effct 
of antiferromagnetism is grossly overestimated. A more 
quantitative approach is to use the renormalized mean- 
field theory by taking into account the constraint a 
little bit more carefully, m^. will be multiplied by a renor- 
malization factor. Instead of pursuing this approach we 
shall use the variational method to calculate the disper- 
sion numerically. Consequently, as shown below, more 
physically reasonable values of parameters are obtained. 
Here, the emphasis is that all the interesting features of 
the dispersion obtained by SCBA is quite consistent with 
the form of Ek- 
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FIG. 1. Energy dispersion of a hole in the t — J model 
obtained by Liu and Manousakis [7] with SCBA on an 20 x 20 
cluster (solid circles), for J = 0.2, and t = 1. Dotted line is the 
fitted result of Ek - Eo with Eo = 5.37, x = 5.47, A = 2.27, 
and rris = 13.2. 

To have a more accurate description of the single hole 
dispersion, we shall use the power method in addition 
to the VMC method. The power method |^^, which is 
essentially the same as the GFMC method, applies many 
powers of Hamiltonian to a trial wave function (TWF) 
to project out the ground state with the same symmetry 
as the TWF. To increase the convergence the method is 
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modified by combining with the Lanczos method, we call 
it the power-Lanczos (PL) method ||2l| , ^ . 

Following the mean- field theory discussed above, we 
can easily construct a variational wave function for a hole 
with momentum k and Sz = 1/2. This function has 
Ne/2 — 1 singlet pairs of electrons and a single unpaired 
electron with momentum k and Sz — 1/2, 

(5) 

where the prime on the summation symbol indicates that 
the momentum k is excluded from the sum if k is within 
the SBZ, otherwise, k-Q is excluded. Notice that is 
essentially the same as C-kxI^'o)- 

The energies obtained by the VMC method for an 
8x8 cluster with J = 0.3 are shown in Fig. 2 as the 
open circles. The variational parameters are A/x — 0.3 
and ms/x = 0.056. Notice that the ground state energy 
for the half-filled lattice (— 22.656J) is subtracted from 
the data. The single-hole state has the lowest energy at 
(7r/2,7r/2). We then apply the PL method to this TWF 
psf to project it onto the lowest energy state. Results 
are shown as the solid circles in Fig. 2. 
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FIG. 2. Energy dispersion of a hole in the t — J model on 
an 8 X 8 cluster, for J — 0.3, and t — 1. Open circles are 
VMC results by using l^*!), solid circles obtained by applying 
the PL method to this wave function. Eik — Eq are plotted 
as the dotted lines. For the lower curve x ~ 6.92, A — 2.71, 
rris = 18.84, Eo = 7 A3, and te// = 0.03 

The results represented by the solid circles in Fig. 2 
are almost identical with that of Fig.l. For both, the 
bandwidth is about 2.2 J, the energy difference between 
states at {tt/2,tt/2) and (tt, 0) is about 0.43 J, and there 
is an EVHS around (tt, 0). However, there is one dif- 
ference. In the SCBA, starting from the classical Neel 
state, the hole only hops between the same sublattice. 
The hole energy begins with / J. States at k and k -I- Q 
have the same energy. Our result shows an energy dif- 
ference between (tt, tt) and (0, 0) and a slight asymmetry 
between energies at k and k + Q. This is due to the 
non-vanishing hopping matrix element between the two 
sublattices which is related with the uniform bond or- 
der X- Kane et al. [p3 have pointed out this difference 
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between a mean-field theory based on a Neel state or a 
RVB state. 

Quantitatively, the difference between Figs.l and 2 
can be accounted for by introducing a coherent nearest- 
neighbor hopping in addition to Ek of Eq. (3), Eik — 
Ek — 2tef f{coskx + cosky). Both the solid and open cir- 
cles in Fig. 2 can be fitted very well by a functional form 
Elk — Eq, where Eq is a constant. Although that solid 
circles represent results much closer to the ground state 
than the results of TWF (open circles), the qualitative 
feature of the dispersion is essentially unchanged. 

Much more accurate numerical results are obtained in 
the exact calculation for a 32-site cluster by Leung and 
Gooding [011 and in GFMC for a 12 x 12 cluster [|. 
In Fig. 3, the 32-site result for J — 0.3 is shown as open 
squares, solid circles are results from GFMC for J = 0.4. 
Both dispersions can be fitted rather well by Eik — Eo 
discussed above. 
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FIG. 3. Energy dispersion of a hole in the t — J model on 
a 32-site cluster (open squares) obtained with exact diagonal- 
ization [10] for J = 0.3, and t = 1. Solid circles are results 
of GFMC [8] for a 12 x 12 cluster with J = 0.4. Dotted lines 
are the fitted Eik — Eq curves. For the lower curve (32 sites), 
X = 5.43, A = 2.43, m, = 17.6, Eq = 7.09, and t,ff = 0.012 
For the upper curve (12 x 12 sites), x = 6.87, A = 2.83, 
nis = 18.4, Eo = 8.9, and t^ff = 0.011. 

The results of our trial wave function shown in Fig. 2 
agree fairly well with the more accurate results of Fig. 3. 
This certainly enhances our belief that the mean-field 
theory or the trial wave function \^i) has captured the 
essential physics of the one-hole state. In particular, we 
notice that the asymmetrical energy dispersion between 
k and k -I- Q are present in both 32 and 144 clusters. 
Although a careful and systematic investigation is nec- 
essary before we can firmly establish this asymmetry for 
larger clusters, this result clearly is in better agreement 
with our mean-field theory than the SCBA calculations. 
As explained above, even for smaU clusters, SCBA §-§] 
results still have the prefect antiferromagnetic symmetry. 

Another interesting result is the observation of the flat 
band region, which is the so called EVHS [^, near (tt, 0) 
in Figs. 2 and 3. Examining Ek or Eik shows that EVHS 
is due to the presence of A or d-RVB order parameter 
and it is independent of the existence of long-range anti- 
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ferromagnetic order. The extra features of the dispersion 
around (tt, 0) depend on the ratio A/x- The d-RVB or- 
der parameter has the largest band gap near (tt, 0). This 
band gap will produce large density of states which is 
reflected by this EVHS. 

In our calculation of the lowest energy state for each 
wave vector k, it is not guaranteed that the state we 
obtained has a finite spectral weight when a hole is pro- 
duced from the half-filled ground state. It is therefore 
necessary to examine the spectral weight defined as 

|(*k|ck,.|$)|', 

where |<I>) is the normalized exact ground state at half 
filling and |^k) is the lowest energy state at momentum 
k in the presence of one hole. Here we use VMC and PL 
method to calculate Zk. The trial functions are |4'o} of 
Eq. (4) and |^i} used to obtain Fig. 2. The results for 
64 sites are plotted in Fig. 4. Open circles are VMC re- 
sults and solid circles are results of first order PL pl| , p3[ . 
We have also included the exact result of 32 sites |lC| ] 
as the open squares. For most wave vectors our results 
have roughly the right magnitude except near k = (0, 0). 
Clearly, we have not yet obtained the exact results for 
64 sites and more higher order power method calcula- 
tions are needed. However, the result is enough to see 
the asymmetry between k and k + Q, which is absent in 
the SCBA. 
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(0.0) (tt.O) (tt,7t) (0,0) (tt.O) (0,7t) 
FIG. 4. Spectral weight Zk as a function of k for the 8x8 
cluster used in Fig. 2. Open circles are VMC results and solid 
circles are by PLl. The exact results of 32 sites are shown as 
open squares. 

So far we have assumed that the hole only has the 
nearest- neighbor hopping matrix element. To explain the 
ARPES !§ resuh, there has been effort (ll]Jl|] to gener- 
alize the model by including the next-nearest-neighbor 
hopping, t' , and next- next-nearest-neighbor hopping, t" . 
Fairly complicated results obtained in these cases can be 
simply understood in our mean-field theory. The basic 
dispersion is determined by the spin interaction. The ad- 
ditional contribution by t, t' and t" is the same as that 
of the ideal gas except with a renormalized magnitude. 
Hence for the t — t' — t" — J model, the one-hole dis- 
persion has the form, E2k = £'ife — At'^^ ^coskxcosky — 
2t"ff{cos{2kx) + cos{2ky)). This is verified in Fig. 5. 




(0,0) (tt.O) (^,tt) (0,0) (m,0) (O,^) 
FIG. 5. Energy dispersion of a hole in the t — t' — t" — J 
model on an 8 X 8 cluster, for J = 0.3, t' — —0.3 and t = 1. 
Circles are for t" — and squares for t" = 0.2. Both results 
are obtained by using first order PL method on |4'i)[23]. Tri- 
angles are the rescaled ARPES data[3]. E2k — Eq are plotted 
as the dotted lines. The parameters are the same as those in 
Fig. 2, except t'^ff = 0.092 for the lower curve, and for the 
upper curve t'^j^f — 0.065 and t'^ff — —0.055. 

In Fig. 5 the first order Lanczos result ||2^ of using 
for the t - t' ~ t" - J model with J = 0.3, and 
t' = — 0.3,<" = is shown as the circles, t' = —0.3 and 
t" — 0.2 as the squares. We have shifted the two sets of 
data so they match at k = (7r/2,7r/2). The lower dotted 
curve shows the fitted result by i?2fc with t'^^j — 0.092 
and t'lfj ~ 0. The squares are fitted with t^^y = 0.065 
and t'l,fj — —0.055. All other fitted parameters are the 
same as those used in Fig. 2. 

Notice that t'^j^ has the opposite sign of t' . This can 
be easily understood in terms of the d-RVB property of 
the TWF. When the hole hops from site ix:iy to site 

+ 1 , + 1 , it rotatcs a nearest-neighbor RVB bond in 
the y-direction to x-direction. Due to the d-symmetry 
of the singlet bond, such a rotation produces a negative 
sign in the hopping matrix element. Similarly for i", a 
nearest-neighbor RVB bond is reversed, and has the 
opposite sign. 

^From Fig. 5 or i?2fe, we can see that the bandwidth, 
which is the energy difference between k = (0, 0) and 
(7r/2, 7r/2), would be reduced if t' is negative and t" — 0. 
But for t" to be positive and of the similar magnitude as 
t', the bandwidth is about the same as the t — J model. 
While for most region of the BZ the effect of t' and t" 
cancels each other, but near (tt, 0) they add up. Hence 
the most significant change between Figs. 2 and 5 is the 
energy increase for the state at (tt, 0). This increase is 
necessary to make the theory in better agreement with 
the ARPES data by Wells et al. (|. The ARPES data are 
the open triangles in Fig. 5. We have rescaled the data 
so that the energies at (0, 0) and (7r2, 7r/2) match that of 
the squares in Fig. 5. Notice that we have not chosen t' 
or t" in order to match the experiments. Even if we had 
used the matched t' or t", the discrepancy between the- 
ory and experiment along (1,0) direction will not be re- 
moved. This can be easily understood by examining E2k- 
The energy at (tt, 0) may be shifted by t' and t", but not 
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at (7r/2,0). Unlike other well known high temperature 
superconductors, Sr2Cu02Cl2 is known to be fairly dif- 
ficult to be doped into a superconductor. We note that 
since (tt, 0) is exactly at the EVHS region with a large 
density of states, some subtle difference between mate- 
rials may become important. The experimental analysis 
is further complicated by the recent observation of the 
d- wave-like gap structure p6| , p7t near (7r,0). More ac- 
curate experimental results for Sr2Cu02Cl2 and other 
high temperature superconductors are needed to clarify 
the discrepancy between theory and experiment. 

In summary, by using a mean-field theory that takes 
into account both staggered magnetization and d-RVB 
singlet we have derived the energy dispersion of a single 
hole in the t — J and t — t' — t" ~ J models. Numeri- 
cal results obtained from exact diagonalization, GFMC 
method and VMC are all in good agreement with the 
dispersion relation. In this theory, occurrence of the 
EVHS near (tt, 0) is due to the presence of d-RVB or- 
der parameter. The dispersion does not have the ex- 
act sublattice symmetry observed in the Neel state in 
which k and k-|-Q are degenerate. This particular fea- 
ture agrees with results obtained by exact diagonaliza- 
tion and GFMC method and disagrees with the SCBA. 
Simple arguments have been provided to account for this 
difference. 

It should be pointed out that there are other quite suc- 
cessful variational studies about the single-hole dis- 
persion. Although the wave functions were constructed in 
different ways, it was recognized that the spin- flip terms 
are essential in obtaining the right physics. The energy 
due to the spin-flip term is exactly what the RVB order 
parameters are designed for. We believe that one of the 
main reasons why so many different calculations men- 
tioned in this paper all have obtained similar results is 
that they all have taken into account the dominant spin- 
flip effect. The RVB theory is the easiest way to take 
into account of this effect right from the beginning. It 
also gives a more intuitively simple interpretation of the 
dispersion and EVHS. 

The EVHS has been observed in many HTS. There 
is not yet a convincing theoretical reason to explain it. 
Although the theory presented here is only valid in the 
presence of a single hole, we can a make a few general 
comments about the underdoped region. In many earlier 
numerical studies the d-RVB state has been known to 
give a good account of the ground state of the t—J model. 
The d-RVB order parameter would produce a gap and a 
large density of state or EVHS at (tt, 0). Such a gap may 
provide a natural explanation for the spin gap and the 
gap observed in ARPES [ p6p^ ] . Theoretical work is now 
in progress to address these issues. 
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